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DRI  Code:  Motivation 


•  Develop  a  fast  in-house  code  for  calculating  the 

Dynamic  Response  Index  (DRI)  injury  metric 

using  test  or  simulation  results  as  input. 

•  Code  should  be  stand-alone  in  nature  and 
should  lend  itself  easily  to  process  automation. 
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Background 


Various  metrics  are  used  to 
predict  the  occupant 
response  and  evaluate  the 
safety  of  vehicle  designs  in 
underbody  blast  events. 

Underbody  blast  events 
cause  a  predominant  risk  of 
thoraco-lumbar  spine  injury. 

The  Dynamic  Response 
Index  (DRI)  has  been  used 
historically  as  a  metric  for 
spinal  compression. 


Vertebral  K 

-  cervical  vertebrae  (7) 

Column 

1 

-  thoracic  vertebrae  (12) 

intervertebral  — S$fc— 
disk  fiK 

-  lumbar  vertebrae  (5) 

sacrum 

coccyx  H _ / 

From  (NATO,  2007). 
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Dynamic  Response  Index  (DRI) 


Measure  of  spinal  injury  risk  that  accounts  for  the  time 
duration  of  a  load. 

Occupant  torso  modeled  as  a  spring-mass-damper  system. 

Calculated  from  maximum  relative  displacement  between 
the  pelvis  and  upper  torso. 


DRI  = 


O^max) 

9 


u>n  =  natural  frequency  (of  spring-mass  system) 
6max  =  maximum  relative  displacement 
g  =  gravitational  acceleration 


Tolerance  level  of  17.7  for  10%  risk  of  AIS  2+  injuries. 
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1-DOF  model 


Takes  anthropomorphic 
test  device  (ATD)  pelvis 
acceleration  or  seat 
acceleration  as  input 
(pelvis  preferred). 

nix  i  (t)  —  FSpring  +  hamper 

mXiit)  =  —kix-L  —  x)  —  c(x!  —  x) 
m(S  +  x)  =  —kS  —  cS 

k  c  . 

6  +  x  = S 5 

m  m 

S  +  x  =  —a>nS  —  2  $o)nS 
x  =  —S  —  2^0)  nS  —  co^S 

-x(t)  =  8  +  2%a>n8  +  <*>^8 


0)0 

 n  max 

S 

•  5  is  the  relative  displacement 
between  the  upper  body  and 
pelvis  =  (xj-x) 

•  L  is  the  damping  coefficient33 
(0.224) 

•  wn  is  the  natural  frequency33 
(52.9  rad/s) 

OX-  k/m 


£ 


C 

2xfmk 
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spinal 

compression 
6  =  x1  -  x 
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1-DOF  Computational  model 


Input  time  and 
acceleration 

•Time  series 
•Triangular  pulse 


Simulation 


Trapezoidal  rule  for 
numerical  integration 
Spring  and  damper  forces 


Acceleration  data  from  physical  test: 
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3-DOF  model 


•  Takes  hull  acceleration  as  input. 

•  Accounts  for  energy  absorption  by 
the  floor  and  seat. 

•  Springs  representing  the  floor  and 
seat  are  piecewise-linear. 

•  The  spring  representing  the  spine  is 
still  linear. 

z3=x3-  x2 
Z2  =  x2-  %x 

Z-i  —  x1  —  X 


m3%3  CO  —  Fspring3  ^dampers 
m3x3(t)  =  -k3(x 3  ~X2)  ~  c3(x3  -  X2) 
m3(z3  +  x2)  =  -k3z3  -  c3z3 
m3z3(t)  =  -m3x2-k3z3  -  c3z3 


Piecewise  Linear  S\ 

2*  Stiffness 

(/) 

LO 

<U 

c 

it 

pring 

to 

bJO 

C 

a. 

Spring  deformation  (m) 

m2%2  CO  —  k’Spring3  Fdamper3  ^spring  2  F damper2 

m2x2(t)  =  k3(x3  -  x2)  +  c3(x3  -  ±2)  -k2{x2  -  x O  -  c2(x2  -  x O 
m2(z2  +  xx)  =  k3z3  +  c3z3  -  k2z2  -  c2z2 
m2z2  CO  =  -m2x1  +  k3z3  +  c3z3  -  k2z2  -  c2z2 


m±Xl  CO  —  ^springs  ^ damper3  ^spring2  ^damper 2 

TTIl^lCO1  =  k2  (x2  —  Xx)  +  c2(x2  —  x-l)  —k1(x1  —  x)  —  c1(x1  —  x) 
m1(z1+x)  =  k2z2  +  c2z2  -  k3z3  -  c1z1 
miZjCO  =  -7 nxx+  k2z2  +  c2z2  _  ^izi  -  ci^i 
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3-DOF  Computational  model 


Input  time  and 
acceleration 

•Time  series 
•Triangular  pulse 


Simulation 


Trapezoidal  rule  for 
numerical  integration 
Spring  and  damper  forces 


Triangular  pulse  input  data  (from 
previously  developed  Excel  code): 


Time  (s)  Time  (s) 
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DRI  Code 


Written  in  Python. 

•  Requires  Python 
2.7+  and  matplotlib 
plotting  library. 

Executed  from 
command  line. 

Allows  several 
optional  arguments. 

Runs  on  Windows, 
Linux,  UNIX,  and  Mac 
OSX. 


Outputs  from  1-DOF 
code  with  triangular 
pulse  input 


Time  (s) 


Pelvis 

—  Upper  body 
Relative 


Time  [sj 
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User  manual 

Explains  input  formatting  and 
output  files  generated. 

Includes  example  command¬ 
line  calls  and  full  test  cases. 

Test  cases  used  to  validate 
code  against: 

—  Previously  developed  Excel  code 

—  Known  DRI  values  for  several 
physical  tests 


a,  •• 

_ /  Home  Insert  Page  Layout  References  Mailings  Review  View 

0  (Body)  -is  a  a  :=  !=  *  *  **  '  I AaBbCcC  AaBbCcC  AaBbCcI  AaBbCcI 

Paste  ,  b  i  ||  ■  *  x,  *'  Ai’lJ!'  A  ■  1  B  •  I  ,A.  •  _  1  Normal  1  No  Spacing  Heading  l  Heading  2  -  chan9' 

“  Styles 

Clipboard  '■*  root  n»  1 1  Paragraph  styles 

3 DOF  Test  Cases 

1)  Input  file:  Test4.txt 

Description:  Triangular  pulse  input  from  previously  developed  Excel  code  for  3DOF 
model  with  time  duration  seconds,  peak  acceleration 

Format: 

DRI  (known): 

(calculated  from  Excel  code) 

Time  Duration  (s)  Peak  Acceleration  (g) 


Command-line  call: 

python  dri3.py  Test4.txt  --skip  4 

Outputs  to  command  prompt: 

Maximum  positive  DRI  value: 

Time  at  which  maximum  positive  DRI  is  reached: 
Maximum  negative  DRI  value: 

Time  at  which  maximum  negative  DRI  is  reached: 
Maximum  lumbar  load: 

Maximum  (absolute  value)  floor  stroke  (m): 
Maximum  (absolute  value)  seat  stroke  (m) : 
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Triangular  pulse  input  data  with  given  time 
duration  and  peak  acceleration: 

DRI  output  from  Excel:  DRI  output  from  Python: 


Acceleration  data  from  physical  test: 

Python  output  DRI  vs.  Time: 


Validation 

•  Validated  against  Excel  code 
for  both  1-DOF  and  3-DOF. 

•  Used  physical  test  results 
with  DRI  calculations 
previously  done  in  other 
software  to  further  validate 
1-DOF  model. 

•  Validated  3-DOF  model 
against  1-DOF  model  by 
setting  very  large  spring 
constants. 
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EARTH  Code:  Motivation 

Error  Assessment  of  Response  Time  Histories  (EARTH) 

-  Compares  time  histories  to  validate  M&S  results. 

-  New  rigorous,  quantitative  tool  for  in-house  VV&A. 

MATLAB  code  delivered  by  the  Automotive  Research 
Center  (ARC)  with  several  papers  but  no  user  manual 

-  (Pan,  2012) 

-  (Sarin,  2008) 

-  (Sarin  et  al.,  2010) 


( 
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Ml  QJ 

Time  -i 

CD 

History  2  § 
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Time  (s) 
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Code 
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MATLAB  R2012b 


EARTH  code 

Combines  existing  measures  and 
algorithms. 

Quantifies  and  separates  error  due  to: 

-  Phase  shift 

-  Magnitude  differences 

-  Topology  (shape)  discrepancy 

Takes  two  time  histories  as  inputs 
along  with  a  few  parameters. 

Outputs: 

-  Plots  of  original,  shifted,  and  warped 
time  histories 

-  Derivatives  of  shifted  and  warped 
time  histories 

-  Error  metrics  for  phase,  magnitude, 
and  topology 

-  Uses  Bayesian  framework  to 
determine  model  confidence  for 
original,  phase-shifted,  warped,  and 
warped  derivative  data. 


PLOTS  APPS 


Hi  &  a 

New  New  Open 
Script  ▼ 

FILE 

XL  CP  New  Variable 

[3  Find  Files  □□  LH 

_/ Open  Variable  ▼ 

Compare  Import  Save  _  i 

Data  Workspace  j Clear  Workspace  ■» 

|  Analyze  Code 

Run  and  Time 

0 f  Clear  Commands  ■< 

CODE 

m  m  i 

►  C:  ►  Users  ►  Dade.Manion  ►  Documents  ►  WA 

►  EARTH  OCT12( 

||  Workspace  f 

Name 

lvalue 

|  Min 

|Max 

n 

±1  CAE 

<4964x1  double> 

3.0150 

3.6050 

id  CAE_temp 

<49604x1  double> 

3.0050 

3.6050 

B  Confidence 

100 

100 

100 

HI  DataType 

<lx4cel> 

B  Ll_magnitude... 

0.0017 

0.0017 

0.0017 

±1  Ll_shape 

0.4511 

0.4511 

0.4511 

U LAMBDA 

205.0184 

205.0184 

205.0184 

±  OutputSelect 

2 

2 

2 

B  PH_0 

0.5000 

0.5000 

0.5000 

±1  SIGMA 

205.0184 

205.0184 

205.0184 

10  Test 

<4964x1  double> 

3.0338 

3.6017 

id  Test  temp 

<49604x1  double> 

3.0251 

3.6188 

Sts 

<49604x1  double> 

25.6100 

31 

Bv 

<49604x1  double> 

3.0050 

3.6050 

cae_der_vvarped 

<5869x1  double> 

-696.4364 

943.4789 

B  cae_eval 

<5869x1  double> 

-696.4364 

943.4789 

i  cae_phase_shif... 

.  <4964x1  double> 

3.0150 

3.6050 

±1  cae_vvarped 

<7710x1  double> 

3.0150 

3.6050 

B  d_eval 

<5869x1  double> 

-172.4711 

178.0700 

B  dt 

1 

1 

1 

E  epsionMethod 

2 

2 

2 

E  epsionhcal 

50.6936 

50.6936 

50.6936 

Warped  Slope  of  Time  Shifted  Histories 
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EARTH  code  validation 


ARC  provided  electrothermal 
battery  model  example: 

•  Test  vs.  simulation  data  for 
terminal  voltage. 

•  EARTH  input  parameters. 

Used  to  ensure  code  was 
working  properly. 

Results  were  consistent  with 
those  of  the  ARC  (Pan,  2012). 


Error  Metric 

Result 

Phase 

1 

Magnitude 

0.0017 

Topological 

0.4511 

Original  Time  Histories 


5000 


Time 


Warped  Time  Shifted  Histories 


0  1000  2000  3000  4000  5000  6000  7000  8000 

Time 
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Acceleration  (m/sA2)  Acceleration  (m/sA2) 


Simple  test  of  EARTH  code 


Original  Time  Histories 


Time  (s) 

Warped  Time  Shifted  Histories 


•  Compared  Excel  output  vs.  Python 
output  for  1-DOF  DRI  model.  Used 
triangular  pulse  input  data  with: 

•  Input  peak  acceleration 

•  Input  time  duration 

•  Down-sampled  from  ~25,000  to 
~1200  data  points  to  reduce 
computation  time. 

•  Very  low  error  across  each  category 
as  expected. 


Error  Metric 

Result 

Phase 

1 

Magnitude 

0.0011 

Topological 

0.0058 

Time  (s) 
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Acceleration  (m/sA2)  Acceleration  (m/sA2) 


Phase  shift  test  of  EARTH  code 


•  Introduced  a  phase  shift  of 
1000  data  points  in  the 
original  vector. 

•  Down-sampled  to  ~800  data 
points,  (phase  shift  of  33) 

•  EARTH  code  recognized  and 
handled  the  phase  shift, 
yielding  low  magnitude  and 
topological  error  again  as 
expected. 


Error  Metric 

Result 

Phase 

34 

Magnitude 

0.0022 

Topological 

0.0103 

Time  (s) 
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Acceleration  (g)  Acceleration  (g) 


Case  1 


Validation  of  EECS  models 


Case  2 


Original  Time  Histories 


Time 


Warped  Time  Shifted  Histories 


Original  Time  Histories 


Warped  Time  Shifted  Histories 


Time 


Time 


Original  data  from 
drop  tower  tests. 


Error  Metric 

Result 

Case  1 

Phase 

3 

Magnitude 

0.4508 

Topological 

0.7635 

Case  2 

Phase 

4 

Magnitude 

0.6779 

Topological 

0.8296 
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Ongoing  work 


Write  EARTH  code  user 
manual  for  in-house  use 
atTARDEC. 

Use  DRI  code  to  help 
evaluate  new  concept 
vehicles  for  DARPA. 

Update  Hybrid  Lumped- 
Finite  Element  (HLF) 
code  (HyperMesh  script 
for  generating  hull 
models)  to  add  3-DOF 
occupant  models. 
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Summary 


1.  DRI 

-  Developed  1-DOF  and  3-DOF  code  in  Python. 

-  Validated  against  Excel  code  and  physical  test  results. 

-  Documented  usage  and  examples. 

2.  EARTH  code 

-  Learned  and  tested  EARTH  code. 

-  Gathered  example  I/O  data. 

-  Applied  to  EECS  Team  data  to  support  in-house 
model  VV&A  efforts. 

-  Documented  code  for  future  VV&A  at  TARDEC. 

3.  Summarized  all  work  in  a  technical  report. 
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